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ABSTRACT 



C^^ , We present the kinematics of a sample of bipolar planetary nebulae (PNe) which 

cover a wide range of observed morphologies and collimation degrees, from bipolar PNe 
with a marked equatorial ring and wide lobes to highly collimated objects. We use 
tJ" ■ an empirical model in order to derive the expansion velocity, collimation degree, and 

r^ ' inclination angle of the PN with respect to the plane of the sky. The equatorial expan- 

^^ . sion velocities measured in the objects in our sample are always in the low to medium 

range (3-16 km s~^), while their polar expansion velocities range from low to very high 
k> , (18-100 km s" ). None of the objects in our sample, even those that show an extreme 

H I collimation degree, seem to be (kinematically) younger than ~ 10^ yr. We compare our 

results with the state-of-the-art theoretical models for the formation of bipolar PNe. We 
find good agreement between the observed expansion velocities and numerical models 
that use magnetic fields with stellar rotation as collimation mechanism. 

Subject headings: planetary nebulae: bipolar — kinematic: planetary nebulae — in- 
dividual(Hen 2-428, K3-58, M2-48, Hen 2-437, K3-46, M3-55, WeSb4, Ml-75, and 
M4-14 
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INTRODUCTION 



The interacting stellar winds model, initially proposed bv lKwok et al.1 (|l978l ). is now the widely 
accepted scenario to explain planetary nebula (PN) formation. However, it has not yet been 
unequivocally determined what is the ultimate mechanism responsible for the collimation of th e 
subclass of PNe that show bipolar or highly collimated morphologies (e.g. iBalick &: Frankl l2002l ) . 
In the context of the interacting stellar winds model, the formation of a bipolar PN requires a 
collimation mechanism. An isotropic fast wind interacting with an aspherical mass loss structure 
formed during the top of asymptotic giant branch (AGB) phase is the most commonly invoked 
scenario, although a collimated fast wind and a spherical AGB structure has also been considered 
(|Lee &: Sahaill2003l ). Despite the fact that models are able to reproduce a wide variety of observed 
morphologies, it is still not clear how the aspherical mass-loss during the AGB (or the axisymmetric 
fast wind) is produced. 

From the numerical standpoint, several variations and additions to the interacting stellar wind 
model have been developed ex tensively in the literature under the name of Generalized Interacting 



Stella r Winds (GISW) model (Balick 



1994 : iMellema fc Frankl Il99,4 llcke et al 



1987 



199 



Icke et al 

1 



198S 



Mellema 



Soker fc Liviolll989l : Frank fc Mellema 



199J). In the GISW model, a fast, tenu- 



ous wind from the central star expands into a slow, dense wind whose geometry is assumed to be 
toroidal. Magnetic fields have also been considered in numerical models; the magnetized wind blown 
bubble model (MWBB) produces an aspherical mass distribution by i ncluding toroidal magneti c 
fields that constrain the outflow and produce jets in the polar direction (JGarcia-Segura et al.lll999l ). 
In ad dition, the precess i on of an episodic j et ha s been used to reproduce po i nt-symmetric morpholo 



gies (ICliffe et al. 



1995 



Steffen fc Loped Il998l : ISoker h Rappaport 



2001 



Garcia-Segura k. Lopez 



2000) • Variations to the magnetic approach to generate the AGB asphe rical density structure in- 



clude using a stellar compariion r apidly rotating around the central star (ICalvet fc Peimbert 



1983 



Garcia- Arredondo &: Frankl l2004l ) . Although models accurately replicate the PNe shapes, the ul- 
timate question still remains: whether the magnetic field, and/or the stellar rotation, required 
to develop an aspherical AGB r nass-loss can be susta ined by a single star or if they require the 
presence of a binary companion (JNordhaus et al.ll2007l ). 



In order to distinguish between the physical processes that may play a role in the process of 
shaping PNe, it is important to have a detailed morphological classifica tion scheme that includes 



morp 
I992I : 



the basic morphological features (jSchwarz et al.lll992l : lGorny et al.lll997l ). and also allows for more 
detailed subclasses (e.g . the presence of multiple shell PNe, multipolar axis) within each group 
( Manchado et al.ll200d ). We should proceed further in the classification and explore the degree of 
collimation observed in bipolar objects, as this most likely reflects the type(s) and strength(s) of 
the physical process(es) involved. Finally, the kinematics of the nebula allows us to recover its 3-D 
structure, as it gives access to the extra-dimension hidden in direct imaging. 

In t his paper, we present h igh resolution echelle long-slit spectroscopy of a sample of 9 PNe 
from the lManchado et al.l (jl996l ) catalog that show highly axisymmetric morphologies with different 



degrees of collimation. Seven of the PNe in our sample are classified as bipolar (K 3-46, K 3-58, 
Hen 2-428, Hen 2-437, M2-48, M3-55, and WeSb4), and the other two (M 1-75 and M4-14) present 
quadrupolar morphologies that are characterized by two pairs of bipolar lobes symmetric with 
respect to two different axes. Two objects in the sample, K3-58 and M4-14, show also point 
symmetric features. In §2 we describe the observations and data reduction, in §3 we outline the 
procedure used in the analysis of the data, and in §4 we present the results of the kinematical fits 
performed to the objects in our sample. Finally, the results are discussed in §5 and the conclusions 
summarized in S6. 



2. OBSERVATIONS 

The images are from " The lAC morphological Catalog of Northern Galactic Planetary Nebu- 



lae" (iManchado et al.lll996l ). The long-slit echelle spectra were obtained in July 1995 and August 
1996 at the 4.2m William Herschel Telescope (WHT) using the Utrecht Echelle Spectrograph (UES). 
A Tektronix CCD detector 1024x1024 pixels, with a spatial scale of 0'.'36 pixel"^ was used. The 
79 line mm~^ echelle was centered on the Ha emission line, covering the spectral range from 6530 
to 6600 A with a dispersion of 0.07 A pixel"^. The slit width was set to I'.'l, providing a spectral 
resolution of 0.14 A which corresponds to 6.5 km s^^. The unvignetted field of view is up to 
160" long. 

The spectra were reduced using standard IRAF tasks for two-dimensional spectra. The images 
were corrected for differential illumination of the slit using sky flat fields obtained with the same 
configuration. The wavelength scale and geometrical distortion were set by a two-dimensional fit 
to a Th-Ar calibration lamp. The wavelength calibration has an accuracy better than 0.004 A (0.2 
km s~^). The sky emission, which includes the geocoronal Ha line, was removed by fitting and 
subtracting the background using a low order polynomial function. 

The observing log is summarized in Table [1] were the common PN name is given in column 
(1), the PN G number in column (2), the date of the observations in column (3), the slit position 
angle (P. A.) in column (4), and the exposure time in column (5). The slit positions go through the 
center of the nebulae, and their P.A. have been chosen to cut across their main symmetry axis. In 
most cases, spectra were obtained at two slit positions along perpendicular directions. 



3. FITTING PROCEDURE 

Bipolar PNe have typically hourglass or butterfly morphology, meaning that they show two 
bipolar lobes connected by a narrow waist. The exact geometry can be very different from one 
bipolar PN to another. These differences can basically be described by the aspect ratio of the bipolar 
lobes and by the relative width of the waist with respect to these bipolar lobes. In addition, the 
inclination with respect to the plane of the sky affects the apparent morphology, while the nebular 



kinematics determines the observed expansion velocity. In order to determine the 3D geometry and 
kinematics of the bipolar PNe in our sample, we have used the empirical model of 



Solf fc Ulrich 



(|l985l ) which was originally developed to interpret the velocity field of the bipolar nebula around 
the symbiotic star R Aquarii. In this model, hereafter referred to as the Solf & Ulrich's model, the 
velocity distribution along a plane going through the main nebular symmetry axis is described by 

v{a) =Ve + {vp-Ve)sin'{\a\) (1) 

where a is the latitudinal angle from the symmetry axis of the nebula, Vg is the equatorial velocity, 
the minimum velocity along the nebular waist, Vp is the polar velocity, the maximum velocity of the 
farthest point from the center of the nebula, and 7 is a geometrical factor that defines the hourglass 
geometry, with 7 ~ 1 for round, bubble- like bipolar lobes, and 7 >> 1 for very elongated bipolar 
lobes. 

One advantage of the Solf & Ulrich's model is the possibility to fit simultaneously images and 
spectra of a given nebula in order to recover the 3D geometry lost by projection effects. As an 
illustration, we show in Fig. [1] simulated images and position- velocity (PV) diagrams of a nebula 
obtained using Eq. [T]for three different inclination angles with respect to the plane of the sky. From 
Fig. [H it is clear that the tilt of the simulated line strongly depends on the inclination angle of the 
object. Therefore, the shape of the PV diagram can be used to derive the inclination angle which 
may only be assessed from the aspect ratio of the central ring under the strong assumption that it 
is circular. Similarly, Fig. [2] shows simulated images and PV diagrams for nebulae with different 7 
factors. In this case, we see that as 7 increases, the shape of the PV diagram becomes narrower and 
the line tilt more apparent. Fig. [2] also shows that it is not possible to determine the inclination 
angle of the nebula from the projected shape of its cylindrical narrow waist for 7 > 9. However, it 
is readily seen that we can determine the inclination from the line tilt in the PV diagram. 

The fitting procedure is an iterative process that is started by modeling simulated images and 
PV diagrams with parameters that are first order estimates. The set of free parameters includes 
the polar and equatorial velocities, Vp and Ue, the 1 kpc kinematical agelj, the inclination angle of 
the PN with respect to the plane of sky, the 7-factor describing the shape of the nebula, and the 
radial velocity. 

For the objects with a clearly visible central ring structure, we can estimate both the equatorial 
expansion velocity and the radial velocity just by measuring the velocity of the two brightness 
maxima at the opposite sides of the ring; the average of the two velocities gives us the radial 
velocity, while the half difference of the velocities give us a lower limit of the equatorial expansion 
velocity. In this case, we can even assess the inclination of the nebula as the axis ratio of the 



^ We use the term 1 kpc kinematical age throughout the text to refer to the age obtained from the kinematical 
fitting assuming that the nebula is located at a distance of 1 kpc. The 1 kpc kinematical age that is determined 
in the model explicitly assumes that the lobes and the waist started to grow simultaneously. The distance- corrected 
kinematical age is the 1 kpc age corrected for the distance to the object. 



equatorial ring is a direct measurement of the inclination angle if we assume that the ring is 
circular. Then, the synthetic nebular image and PV diagram are plotted over the Ha or Hq;+[N ll] 
image and [N ii] spectr4j, respectively. The initial estimates of the parameters are changed until 
we reach an optimum fit when the synthetic PV diagram passes through the intensity maxima in 
the spectral data and the simulated image traces tightly the nebular geometry. 

We note that the geometry of the nebula can be reproduced with different sets of values for the 
free parameters. However, the degeneracy is broken by fitting the geometry of the object together 
with the spectral line shape. The accuracy of this iterative process is hard to quantify and it has 
to be determined for each object individually. 

The parameters of the best fits are listed in Table 2 and Figs. [3] to [8] show the images and 
PV diagrams overplotted with our best fits to the data. The velocities in the horizontal axis of 
Figs. [3] to [8] have been converted to the LSR system. The images are shown at the same scale as 
the spectra and have been rotated to allow a direct comparison of the kinematical features along 
the spatial direction. 



4. RESULTS: THE FITS TO THE IMAGES AND SPECTRA 

In the following we provide a detailed description of both the data and the best fit parameters 
obtained for each object. 

Hen 2-428 

Hen 2-428 (Fig. [3|) is a bipolar PN with a noticeable equatorial ring and two open hourglass 
bipolar lobes of which the Northern one is brighter and m ore extended. The cen tral star, unusually 
bright for bipolar PNe, has a known binary companion ( Rodriguez et al.ll200ll ). For this nebula. 



we acquired long-slit echelle spectra along PA's, 77° and 157° (see Fig. [3]). We note that the 
spectrum taken along P.A. 77° is off-center and passes through the bright knots in the equatorial 
ring. Therefore, the basic spectral fitting relies on the spectrum at P.A. =157° which has adequate 
signal-to- noise ratio in the central region and the brightest Northern lobe to constrain the fit. From 
our best fit, we get an equatorial expansion velocity, Ve, of 16 km s~^ , and a polar velocity, fp, of 



80 km s . A value of fe=15 km s was obtained for this object bv lRodriguez et al.l (|200ll ). 



Hen 2-437 

Hen 2-437 (see FigH] [Top]) is a very elongated PN that does not show a central ring structure. These 
morphological features make very difficult to constrain its equatorial velocity and to reproduce the 
spectral shape with the Solf &: Ulrich's model. Since there is almost no tilt in the central maxima 
of the [N ii] echellogram along PA 77°, we assume that the nebula has a small inclination angle, 
lying almost in the plane of the sky. 



For the spectral fits, we use the [N ll] emission Une because of its relative brightness and small thermal broadening. 



The geometry of the nebula is thus constrained by its morphology, and the nebular contours 
can be reproduced relatively well by using an extremely high 7 factor. Given the difficulties in 
reproducing the spectral shape of Hen 2-437, it was not possible to derive firm results. We give a 
set of the most likely values: the inclination angle is small, within 2°, the polar velocity can be 
constrained to be between 50 and 100 km s~^, and the equatorial velocity has to be small, of the 
order or lower than 10 km s~^. 

K3-46 

K 3-46 (Fig. S] [Bottom]) has a well defined hourglass geometry with a prominent equatorial ring 
and a marked waist. For this nebula, we obtained very low expansion velocities from the fit to the 
PV diagram, with v^, < 3 km s~^and Vp ~ 18 km s~^. We note that it was not possible to find 
an optimal fit, and Fig. [4] [Bottom] indeed suggests that our simulated fit would require a larger 
central ring. We also note that the kinematical age of K 3-46 is of 24,000 yr which is, as we will 
discuss later, the highest value obtained in our sample. 

K3-58 

At a first glance, K3-58 shows a classical bipolar morphology (Figj5]), but its conspicuous ring 
is rather irregular and the bipolar lobes show point-symmetry. All these components are clearly 
detected in the [N 11] echellogram obtained along P. A. 90°, with the brightest emission peaks located 
on the equatorial ring. Inside the main lobes, a secondary structure reaches approximately up to 
one half of the length of the main lobes. This secondary structure, that is clearly present in the 
spectral data as bright emission peaks, may have originated from a later ejection to the one that 
shaped the main lobes. From the fit we obtain t;p=38 km s~^ and Ve=12 km s~^. Interestingly, the 
ratio between the equatorial and the polar velocities is relatively small for this object. 

Ml-75 



M 1-75 (Fig. [6]) is the first object classified as quadrupolar (JManchado et al.lll996l ). Ml-75 with 
its complex morphology and extremely high He abundances and N/0 ratio is a candidate for a 
peculiar stellar evolutionary path (Guerrero et al. 1996). The nebula does not show a circular 
equatorial ring, but an irregular annular structure at the central regions. Since the equatorial ring 
is clearly not circular, it cannot be used to assess the inclination angle of the nebula with respect 
to the plane of the sky. 

The spectra for M 1-75 was obtained along the main symmetry axis of the two pairs of bipolar 
lobes, namely along P. A. 150° and 176°. The former spectrum (Fig. [6] [Top]) proved more useful 
to guide our spectral fit, as it provides information for both pairs of bipolar lobes. The data taken 
at the second position angle has been used for control purposes only and is shown at the bottom 
of Fig. El 

The large lobes, along P. A. 150°, can be well fitted by using 7=6 and an inclination angle of 
87° as determined from the spectra. Since these lobes seem to lie close to the plane of the sky, a 
change in just 2° is very noticeable for such an elongated shape, and the inclination angle can be 
constrained well. The expansion velocities obtained for the largest pair of lobes are relatively low. 



with Ue=8 km s ^ and fp=55 km s ^. 

The kinematics of the smah lobes cannot be constrained as accurately because their geometry 
in the image is not well defined. We have assumed, for simplicity, that the small and large pairs of 
lobes have identical morphologies, and therefore we use the 7 factor determined for the large lobes, 
7=6. The range of inclination angles we then obtain is between 60° and 70°. Using an inclination 
angle of 65°, we obtain expansion velocities of t'p=45 km s~^ and Ve=12 km s~^. 

M2-48 

M2-48 (see Fig. [7] [Top]) does not show an equatorial ring, but a pinched waist. The pair of 



collimated lobes ends in a pair of bow-shocks that were analyzed in detail by iLopez-Martin et al 



(|2002l ) . The collimated morphology of the bipolar lobes of M 2-48 requires to use a high 7 factor 
(~ 8). The narrow central waist of M2-48 does not allow to determine the inclination angle from 
the geometry of the ring and, therefore, as mentioned in §3, we have to rely on the line tilt in the 
observed PV diagram. We reproduce well the geometry of both image and spectra by using an 
inclination angle of ~ 80°, with the North-Eastern lobe receding from us. We derive a radial velocity 
of 16 km s~^, a Ve=W km s~^, and a Vp of 100 km s~^ which is the largest value in our sample. 
The inclination ang l e and expansion velocities obtained are consistent with the determination by 



Lopez-Martin et al.l (120021 ). 



M3-55 

Despite being very faint, M 3-55 displays a clear symmetric shape (Fig. [7] [Bottom]). We obtain for 
M3-55 a low expansion velocity, lower than 6 km s~^ in the equator and about 19.5 km s~^ along 
the polar directions. M3-55 also has the lowest 7-factor of all the objects within our sample which 
is not surprising given that it has a relatively wide ring and round lobes. 

M4-14 



Manchado et al.l (|l996l ) classified M4-14 (Fig. [8] [Top]) as a PN with quadrupolar morphology 
(based on imaging). Here, we also remark the noticeable point-symmetry of its bipolar lobes. To 
fit the cylindrical waist of M4-14, a high 7 factor (7=5) is needed. From this fit we obtain v^ = 11 



km s ^ and Vp =65 km s ^. 



M4-14 has a [N 11] /Ha ratio among the highest found in PNe, indicative o f a high nitrogen en- 



hanc ement and high N/0 ratio, that are known to be correlated with bipolarity (jPeimbert &: Torres-Peimbert 



19831 ). The chemical abundances, expansion velocity, and geometry factor 7 make M 1-14 and M 1- 



75 very similar. 

WeSb4 

WeSb4 is a large, somehow diluted object, that shows one of the most irregular morphologies 
among the objects in the sample (Fig. [8] [Bottom]). The narrow-band image does not display a 
clear hourglass shape that, on the other hand, is evident in the much deeper spectrum. In the 
spectrum we detected weak, and more extended emission than in the optical images. From the 
image and [N 11] echellograms along P. A. 69° and 159°, we derive an inclination angle of 50°, and 
a Ve of 14 km s~^, and Vp of 95 km s^^WeSb4 has one the highest Vp in our sample. 



5. DISCUSSION 

We summarize in Table [2] the parameters obtained from the best fits to the data. Column (2) 
gives the size of the lobes and waist as measured from the maximum extension of the 3a contour 
levels extracted from the images obtained by Manchado et al. (1996). The polar and equatorial 
expansion velocities are given in column (3) and (4), respectively, the 1 kpc kinematical age is 
given in column (5), the inclination angle in column (6), and the 7 factor in column (7). Below 
the values obtained from our best fit (when the synthetic PV diagram passes through the intensity 
maxima in the spectral data and the simulated image traces tightly the nebular geometry), we have 
listed the range of values that still provide reasonable matches (when we can still provide a fit that 
passes through the main features) to the images and spectra in order to provide an estimate of the 
uncertainties in our fits. 

The equatorial velocities range from very low values (3 km s^^ for K 3-46) to typical expansion 
velocities (16 km s^^ for Hen 2-428). None of the objec ts in our sample has high , ~40 km s^^. 



equatorial expansion velocities such as the ones found by lCorradi &: Schwarzj (|l993l ) in the bipolar 
PNe CTS 1 and Hen 2-84. The polar velocities cover the whole range, from very low (18 km s~^ in 
K3-46) to medium/high (100 km s"^ in M2-48). 

The PNe analyzed in this paper can be classified into three different groups, according to 
their morphology and kinematics. The first group, formed by Hen 2-428, K3-46, K3-58, M3-55, 
and M4-14 consists of bipolar PNe with notable central rings. For these objects, we have a direct 
estimate of their equatorial expansion velocities measured from the spectral emission of the ring (as 
it was described in §3). A second group is formed by the highly collimated objects Hen 2-437 and 
M2-48 having the highest polar expansions and elongated bipolar morphology typical mostly for 
younger PNe. Finally, the third group is formed by the somewhat deteriorated objects with poorly 
defined morphologies; M 1-75 and WeSb 4. The PNe belonging to this third group are expected to 
be more evolved since they have the highest kinematical ages. Although K3-46 also has a large 
kinematical age, we did not include it in this group of evolved objects. We suspect that its large 
kinematical age is probably a consequence of deceleration in the course of its evolution. 



5.1. PN Ages 

Our spatio-kinematical study provides a direct estimate of the kinematical age that can be 
used to assess the nebular age. It is important to keep in mind in this comparison the very likely 
possibility of acceleration or deceleration of the nebular material due to the complex interaction 
between the ionization and dynamics of the shell driven by the hot bubble. A s a result, kineniatica l 



ages often do not match the age of the central stars, as shown theoretically bv lVillaver et al.l (|2002l ). 
who proved that kinematical ages overestimate the age in young nebulae and underestimate it for 
evolved ones. Moreover, the definition of PN age is a rather tricky concept. It can be considered 
that a PN is born when the central star supplies enough photons capable of ionizing the nebula. 



However, the kinematical age is a dynamical concept that tell us when the gas started moving. At 
the time of nebular ionization the gas is already moving and the fotoionization its elf is expected 



to change the gas dynamics. In addition, the gas velocity as inferred from models IVillaver et al 



(|2002. ) is not constant. If the gas has been accelerated the kinematical ages are underestimated 
and the opposite is true if the gas has suffered deceleration. The age obtained in Table 2 is one of 
the parameters of the fitting and it is obtained under very simplistic assumptions, namely that the 
lobes and the waist were formed at the same time and they have been moving at a constant velocity 
since then. Therefore, the ages determined this way are just an order of magnitude approximation 
to the time since the formation of the nebula (both dynamically and from ionization). 

The kinematical model assumes the same age for the lobes and the waist of a nebula. As a 
consequence the larger lobes always have larger expansion velocities. Note that under this simplistic 
assumption we are excluding the possibility that the lobes were formed before the waist. 

We provide in Table [2] the 1 kpc kinematical ages. Here, we have estimated the "real" ages 
by multiplying the 1 kpc kinematical ages derived from our fits by the individual distances. It is 
well known that PNe distances are poorly determined, however they are a necessary parameter 
for the age estimation. It is important to note that the ages determined in Table 3 are the ones 
given in Table 2 but scaled to the distance to the nebula. We have used the distances obtained 
from the Galactic rotation curves of Burton (1974) when possible, as these distances are considered 
to be the most reliable. In some cases, our data was out of the range covered by the rotational 
curves, and therefore the distance estimate was not viable through this method. This also applies 
for M2-48, for which this method gives ambiguous results, and for Ml-75, for which the distance 
derived from this method (5.4 Kpc) is unreasonable and gives a too long kinematic age for a PN. 
When it was not possible to estimate the distance from the rotational curves, we have used the 
stati stical distances from Cahn et al (1991) when available. Otherwise we have used the distance 



from ICahn &: Kaleij (|l97lh or, as a last resort, an average of available distances in the literature. 



We note that the values given bv iMaciell (| 19841 ) are systematically lower than the values given by 



others. Therefore, Maciel's distances have not been used for the averages. 

In Table[3l column (1) gives the PN name, column (2) the radial velocity, column (3) the kine- 
matic age at 1 kpc obtained from the model fitting, column (4) the statistical distances taken from 



Acker et al.l (jl992l ). column (5) the distances estimated from the galactic rotation curve (Burton 
1974), and column (6) the distance used to estimate the kinematical age which is given in column 

(7). 

The distance-corrected kinematical ages range between ~5,000 and ~20,000yr. Intuitively, it 
is expected that younger PNe would have better defined morphologies than the older ones. This 
appears to be the case for almost all the objects analyzed. Hen 2-428, M2-48, M3-55, and M4-14 
are relatively young, with kinematical ages ~5,000 yr, and all have well defined morphologies. 
M 1-75 and especially WeSb4 are older, with kinematical ages ~10,000 yr, and their morphologies 
are not sharp. 



10 



In old objects such as WeSb 4, ionization-driven instabilities which act on a timescale com- 
parable to the kinematica l age of the nebula might be responsible for the development of irregular 
shapes (jciiffe et alJll995l ). 



There are two objects in our sample for which a relation between their age and their morphology 
is not straightforward. K3-46, the oldest PN in our sample (~20,000 yr), and K3-58, a relatively 
old PN, have sharp morphologies. One possibility is that the distances used are wrong, which would 
not be a surprise given the high uncertainties involved in PN distance determinations. Another 
possibility for K3-46 is that it has experienced deceleration, in which case the age given in the 
Table 3 would be an underestimation of the real age of the object. 

Adopting a distance of 1 kpc to Hen 2-437, this is the youngest PN in our sample, with a 
kinematical age be tween 750 and 2000 yr. This age is consistent with those of high collimated PNe 
(|Lee fc Sahailbood ^. 



5.2. Comparison with Numerical Models 

Quite often the comparison between theoretical models and observed objects is based only on 
a morphological match since this the only available information in most cases. From our experience 
in matching published models to data, we have found that it is possible to find models that can 
simulate the morphology of a PN but with velocities different to those measured; in this case we 
need to know whether the model can simultaneously reproduce the morphology and kinematical 
properties that, in the case of bipolar PNe, are parameterized by the expansion velocity, the physical 
nebular size, and geometry (the 7 factor). Note that, in the model of Solf & Ulrich, the expansion 
of a bipolar PN is homologous, i.e., in its expansion, the nebula keeps its proportions, as well as 
the ratio between the polar and equatorial velocities. 



Only the models by iGarcia-Segura et al.l (jl999l ) and iGarcia-Segura &: Loped (J2000l ) provide 
enough information on the evol ution of the morphol o gy for different collimation parameters for us 
to compare to our data. In the iGarcia-Segura et al.l (119991 ) models (hereafter, GS99), the velocity 
ratio between the equator and the poles depends on the shaping mechanism, while the absolute 
values of the velocities depend on the initial velocity of the slow wind. Therefore, when evaluating 
the ability of the theoretical models to match our observations, we are more interested in comparing 
the nebular shapes and velocity ratios rather than the absolute velocity values. 

To compare the models with the observations, we need a parameter that can quantify the degree 
of collimation. The 7 factor in the Solf & Ulrich geometrical model is an appropriate parameter to 
account for the degree of collimation observed. The shape in the GS99 models is controled using 
different values of a (the ratio of the magnetic to the kinetic energy density in the fast wind) and 
(ratio of the stellar rotational velocity to the critical breakup velocity). 



In order to allow a better comparison between the models and the parameters derived from 
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the observations we have determined the Solf & Ulrich's parameters (7 factor) that correspong to 
the GS99 simulated cases. We find that low a and 17 values produce shapes that are similar to 
those generated by low values of 7 (7 < 1) (see Fig. 5 in GS99). The GS99 cases with low a and 
high Q values are equivalent to medium 7 values, while high a and high 0, values result in very 
collimated objects, with high 7 values. 

For categories J~K in GS99, the morphology develops using only rapid star rotation rates, 
while for categories Q-V the additional help of magnetic fields is required. Hen 2-428, K3-46, and 
M3-55, which have the lowest 7 (1, 0.9, 0.6, respectively) factors in our group, compare well to the 
categories J-K of GS99. Both the observed objects and the corresponding theoretical simulations 
are reproduced with 7 factors around 1 and, although our PNe do not have the same absolute 
expansion velocities as the GS99 cases, this is not a critical issue as lower expansion velocities 
in GS99 magneto-hydrodynamic (MHD) simulations could be obtained by choosing lower initial 
velocities for the slow wind. In fact, the ratio of expansion velocities for K3-46 and Hen 2-428 
situates them in categories Q-R, but their 7 factor and morphology place them into the J-K 
group. It must be noticed that the central star of Hen 2-428 is in a binary system (Rodriguez et al. 
2001). It is tantalizing to consider that the presence of the binary companion may have influenced 
the nebular shaping, resulting in a bipolar PN with enhanced polar velocities with respect to the 
bipolar PN that would have been collimated by a single, rotating star. 

In general, the expansion velocities in GS99 are rather high compared to those observed in our 
sample. However, M 1-75, M4-14, and WeSb4 have similar expansion velocity ratios (i.e. equatorial 
to polar expansion velocity) and morphologies as the PNe in categories Q~R in GS99. These 3 PNe 
have 7 factors between 5 and 6. In fact, M4-14 even coincides in its kinematical age and value of 
its expansion velocities with the example given by GS99. 

M2-48 is more collimated than the previous objects and fits categories S-T, while Hen 2-437 
would go even higher, to the class U in GS99. Such highly collimated objects, with lobes of wedge- 
shaped polar regions, are not well reproduced by the Solf & Ulrich's model and therefore our derived 
velocities carry a large uncertainty. 

Careful inspection of the images of M 2-48 and Hen 2-437, which have the highest collimation 
factors in our sample, reveals traces of a dusty disk in the equatorial region, a feature not visible 
in the I0W-7 objects. Similar features were already detected in various young, highly collimated 



objects (e.g. ISahai et al.l 120061 ) and suggest the existence of thick equatorial disks. 



Finally, we are left with K 3-58 a rather peculiar object since it has an expansion ratio and 7 
factor that lie outside of the cases modeled by GS99. K 3-58 presents a cylindrical equatorial ring 
which is wider than the one obtained for the J-L categories of GS99, and it shows lobes that are 
rounder than those in the R-V classes. The best fi t to th e morphology and kinematics of this object 



can be found in Fig. 5 of iGarcia-Segura &: Loped (J200d ). The nebula simulated there has almost 



identical expansion velocities to those we measure. In this model, the point-symmetric shape is 
accomplished with a tilt between the magnetic collimation axis and the bipolar wind outflow. 
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To summarize, we find that objects with low 7 values generally compare well to the GS99 
models which are mainly shaped by stellar rotation. The possible exceptions are PNe formed in 
binary systems that will present low 7 but somewhat higher expansion velocities (Hen 2-428), while 
high 7 objects show high polar expansions and agree well with the simulations obtained using 
strong magnetic fields and high rotational velocities. 



6. CONCLUSIONS 

Although all of our objects are bipolar/quadrupolar, they differ significantly in their mor- 
phology and kinematics. Morphologies range from bubble-lobed (Hen 2-428) to highly collimated 
(Hen 2-437) PNe. The sharpness of the nebular shapes also varies from well defined objects to these 
with somewhat deteriorated shapes (WeSb4). This variety is reflected in the range of geometrical 
7 factors (from 0.6 to 20) which might differentiate between the physical processes that originate 
them. 

We find that the objects from our sample present as well a variety of expansion velocities, 
from low (3 km s~^) to medium (16 km s~^) equatorial expansions, and from low (18 km s~^) to 
medium/high (100 km s~^) polar expansion velocities. The disagreement between the spectral data 
of K 3-46, hardly revealing signs of equatorial expansion, and its image, showing a wide central ring 
indicate the possibility of deceleration in this nebula. This deceleration may reveal the interaction 
of the nebular material with a dense equatorial disk. 

The estimates of kinematical ages, derived using distances inferred from the Galactic rotation 
curve or otherwise statistical distances, vary from middle age to old with possible significant errors 
originated from distance errors and non-uniform expansions. The data agree rather well with the 
state-of-the-art theory for PN collimation, however, we cannot exclude the origin of the possible 
shaping mechanisms i.e. whether rapid star rotation, and/or magnetic fields are originated by 
single or binary systems. We suggest that the 7 factor used to fit the Solf & Ulrich's model could 
roughly indicate which shaping process is actually at work. 
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Fig. 1. — Synthetic images and position- velocity (PV) diagrams for different inclination angles for 
a bipolar PN with 7=1. 



Fig. 2. — Synthetic images and position-velocity (PV) diagrams for a bipolar PN at an inclination 
angle of 30° with different values of 7. 



Fig. 3. — Kinematical fitting of the Hen 2-428 for a slit position angle of 167° (top) and fo r a sli t 
position angle of 77° (bottom). The narrow-band image, extracted from the lManchado et al.l (119961 ) 



catalogue, has been rotated to the position angle of the slits along which the spectra were obtained. 
The angular scale of images and spectra have been matched to allow an easy comparison of images 
and spectra. Also note that the PV diagrams have been corrected for the local standard of rest 
(LSR) velocity. In order to display the complete structure of the spectral line, the contrast of the 
spectra is adjusted, so the faintest parts can be seen. 



Fig. 4. — Same as Fig. 3 for Hen 2-437 for a slit position angle of 77° (top) and K3-46 for a slit 
position angle of 106° (bottom). In the later, the central part of the spectral line is shown saturated, 
and white contours have been used to show levels of the same intensity. 



Fig. 5. — Same as Fig. 3 for K3-58 for a slit position angle of 90° (top) and for a slit position angle 
of 14° (bottom). As for K3-46, white contours are used for the brightest regions of the spectral 
line. 



Fig. 6. — Same as Fig. 3 for M 1-75 for a slit position angle of 150° (top) and for a slit position 
angle of 176° (bottom,). As for K3-46, white contours are used for the brightest regions of the 
spectral line. 



Fig. 7. — Same as Fig. 3 for M2-48 for a slit position angle of 66° (top) and M3-55 for a slit 
position angle of 57° (bottom,). As for K3-46, white contours are used for the brightest regions of 
the spectral line. 



Fig. 8. — Same as Fig. 3 for M4-14 for a slit position angle of 40° (top) and WeSb4 for a slit 
position angle of 159° (bottom). As for K3-46, white contours are used for the brightest regions of 
the spectral line of M4-14. 
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Table 1. Spectroscopic Observations 



PN NAME 


PN G Number 


Date 


SlitP.A. 

(°) 


ExposurcTime 


Hen 2-428 


49.41+2.48 


1995 Jul 15 


77 
167 


1800 
1800 


Hen 2-437 


61.40+3.62 


1995 Jul 15 


77 
167 


1800 
1800 


K3-46 


69.21+3.81 


1995 Jul 16 


16 
106 


1800 
1800 


K3-58 


69.68-3.94 


1996 Aug 07 


14 
90 


900 
1800 


M 1-75 


68.86-0.04 


1996 Aug 07 


150 
176 


1800 
1800 


M 2-48 


62.49-0.27 


1996 Aug 07 


66 


1800 


M3-55 


21.74-0.67 


1995 Jul 15 


138 
57 


1800 
1800 


M 4-14 


43.09-3.04 


1996 Aug 07 


40 
87 
153 


1800 
1800 
600 


WeSb 4 


31.91-0.31 


1995 Jul 14 


69 
69 
159 


1800 
1800 
1800 
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Table 2. Results of Spatio-Kinematical Modeling 



PN Name 


Ha size^ 


Vp 


Ve 


Kin. Age 


Incl. Angle 


7 




["] 


[kms-i] 


[kms-1] 


[yr] 


[°] 




Hen 2-428 


63x18 


80 


16 


2400 


-75 


1 






[73,85] 


[14,17] 


[2400,2500] 


[-74, -76] 


[0.9,1.1] 


Hen 2-437 


45x4.6 




(5) 




-89 


20 






[50,100] 


< 10 


[750,2000] 


[-88, -90] 


[10,45] 


K3-46 


81x36 


18 


3 


11000 


-70 


0.9 






[22,32] 


[1,4] 


[7000,12000] 


[69,72] 


[0.6,1.2] 


K3-58 


23.0x12.7 


38 


12 


1800 


-60(-57) 


2.5 (2.2) 






[32,42] 


[10,15] 


[1500,2200] 


[-55,-63)] 


[2.1,2.7] 


M 1-75 


69x232 


55 


8 


2700 


87 


6 






[45,65] 


[6,10] 


[2300,3400] 


[86,88] 


[5.5,6.5] 


M 1-75^ 




45 


12 


2400 


65 


5.5 






[25,50] 


[6,13] 


[2000,5200] 


[60,70] 


[5.3,5.7] 


M2-48 


42.6x5.7'* 


100 


10 


1160 


-79 


8 






[75,120] 


[8,11] 


[800,1600] 


[-73, -80] 


[7,9] 


M3-55 


11.1x8.2 


19.5 


6 


1800 


40 


0.6 






[17.5,23.5] 


[4,8] 


[1000,2500] 


[30,50] 


[0.4,0.9] 


M4-14 


27x8.5 


65 


11 


1500 


38 


5 






[44,76] 


[7,11] 


[1400,2500] 


[37,41] 


[4.7,5.2] 


WeSb4 


77x18"* 


95 


14 


3400 


50 


6 






[80,110] 


[13,17] 


[2750,3700] 


[48,53] 


[5.8,6.2] 



* The size is of the lobes and waist 
2Ring size 24.0x13.1 
^For the smaller pair of lobes 
^Size in the [N II] image 
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Table 3. Estimated Kinematical Ages 



PN Name 



( kins^-*-) 



1 kpc Kinematic 
Age 

(yr) 



statistical 

Distance 

(Kpc) 



Reference 



Rotation 

Curve Dist. 

(Kpc) 



Adopted 

Distance 

(Kpc) 



Estimated 
Age 

(yr) 



en 2-428 


72 


3-46 


66 


3-58 


21 



Ml-75 



M2-. 



M3-55 



M4-14 



WeSb4 



16 



30 



49 



2400 

9000 
1800 

2700 



1160 



1800 



1500 



3400 



2.7 


CaKa71 


- 


1.7 


Ma84 




2.15 


CaKa71 


- 


6.6 


CaKa71 


- 


5.7 


Ma84 




2.6-3.7 


CaKa71 


5.3 


3.9 


Ca76 




2.8 


Ac78 




3.21 


Da81 




2.40 


AGNR84 




7.0 


Ma84 




3.89 


CKs91 




3.42 


CaKa71 




1.6 


Ma84 


1.5-7.7 


6.97 


CKs91 




3.56 


CaKaVl 


2.8 


1.4 


Ma84 




3.13 


Ac78 




3.7 


CaKa71 


3.7 


1.6 


Ma84 




2.97 


Ac78 




6.69 


CKs91 





2.2 

2.2 
6.2 

3.4 



4.7 



4.2 



2.8 



3.7 



4.7 



5280 

19350 
11070 

9099 



4872 



5040 



5550 



15980 



Table 3. References: [CaKa71] - Cahn & Kaler (1971); [Ma84] - Maciel(1984); [Ca76] - Cahn (1976); 
[Ac78] - Acker (1978); [Da82] - Daub (1982); [AGNR84] - Amnuel et al. (1984); [CKs91] - Cahn et al 

(1991) 
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